Diffusion-Limited Coalescence, A + A # A, with a Trap 
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We study diffusion-limited coalescence, A+A ^ A, in one dimension, and derive an exact solution 
for the steady state in the presence of a trap. Without the trap, the system arrives at an equilibrium 
state which satisfies detailed balance, and can therefore be analyzed by classical equilibrium methods. 
The trap introduces an irreversible element, and the stationary state is no longer an equilibrium 
state. The exact solution is compared to that of a reaction-diffusion equation — the habitual 
approximation method of choice. The exact solution allows us to determine the rate coefficients in 
the reaction-diffusion equation, without appealing to renormalization group techniques. 
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I. INTRODUCTION 

Non-equilibrium kinetics of diffusion-limited reactions 
has been the subject of much recent interest |j]-Q|. In 
contrast to equilibrium systems — which are best an- 
alyzed with standard thermodynamics — or reaction- 
limited processes — whose kinetics is well described by 
classical rate equations j^U — there is no general ap- 
proach to non-equilibrium, diffusion- limited reactions. 

In this letter we study a diffusion-limited coalescence 
process in one dimension: A + A ^ A, which can be an- 
alyzed exactly ]j|JlO[|. When a trap is introduced, the re- 
sulting steady state is a non-equilibrium state. We derive 
an exact description of this state and compare it to the 
prediction from a reaction-diffusion equation — the stan- 
dard approximation method of choice. The exact solution 
allows us to determine the appropriate rate constants of 
the reaction-diffusion equation directly. Normally, this 
feat requires a renormalization group analysis. 

The rest of this letter is organized as follows. In Sec- 
tion [n] we present a lattice model of diffusion-limited 
reversible coalescence, along with the exact method of 
analysis; the method of Empty Intervals, also known as 
the method of Inter-Particle Distribution Functions. The 
stationary state in the presence of a trap is derived and 
analyzed in Section III . In Section [i^ we compare the ex- 
act solution to that of a reaction-diffusion equation, and 
devise strategies to determine the appropriate rate coeffi- 
cients. We conclude with a discussion and open questions 
in Section [v|. 



II. REVERSIBLE COALESCENCE 



Our model ||,|10| is defined on a one-dimensional lattice 
of lattice spacing a. Each site is in one of two states: oc- 
cupied by a particle A, or empty. Particles hop randomly 



into nearest neighbor sites, at rate D/a 2 . A particle may 
give birth to an additional particle, into a nearest neigh- 
bor site, at rate v/2a (on either side of the particle). If 
hopping or birth occurs into a site which is already oc- 
cupied, the target site remains occupied. The last rule 
means that coalescence, A + A — ► A, takes place im- 
mediately upon encounter of any two particles. Thus, 
together with hopping and birth, the system models the 
diffusion-limited reaction process A + A ^ A. 

An exact treatment of the problem is possible through 
the method of Inter-Particle Distribution Functions 
(IPDF). The key concept is E ni7n (t) — the probability 
that sites n, n + 1, • • • , m are empty at time t. The prob- 
ability that site n is occupied is 



Prob(site n is occupied) = 1 — E n ^ n 



(1) 



The event that sites n through m are empty (prob. E n ^ m ) 
consists of two cases: site m + 1 is also empty (prob. 
E n ,m+i)i or it is occupied. Thus the probability that 
sites n through m are empty, but site m + 1 is occupied 
is E n<m — E n _ rn+ i. With this (and with a similar rule for 
when the particle is to the left of the empty segment) one 
can write down a rate equation for the evolution of the 
empty interval probabilities: 
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Eq. (||) is valid for m > n. The special case of m = n 
yields the boundary condition 
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(3) 



The fact that the {E n ^ m } represent probabilities implies 
the additional condition that E n ^ m > 0. Finally, if the 
system is not empty then E n _ m — > as n — > — oo and 
m — > oo. 

In many applications, it is simpler to pass to the con- 
tinuum limit. We write x = na and y — ma, and replace 
E n , m (t) with E(x,y;t). Letting a — > 0, Eq. (0) becomes 
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with the boundary conditions, 



E{x, x; t) = 1 , 
£(z,y;t) > , 
lim y; i) = . 



(5) 
(6) 
(7) 



The concentration of particles is obtained using 
Eqs. (|l|) and (^), and passing to the continuum limit: 



(8) 



It can also be shown ]g|Jl0| that the conditional joint 
probability for having particles at x and y but none in 
between, is 



P(x,y;t) 



d 2 E(x,y;t) 
dx dy 



(9) 



Given a particle at x, the probability that the next near- 
est particle is at y, i.e., the IPDF, is p(x,y;t) = (l/c)P. 



III. STEADY STATE WITH A TRAP 

The steady state of Eq. (Q), with the boundary condi- 
tions (|) - (|), when dE/dt = 0, is 



E eq = e-^v-i 



(10) 



This corresponds to the equilibrium concentration of par- 
ticles (using Eq. |h 
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(11) 



Another trivial solution is E(x, y) = 1: it represents a to- 
tally empty system. But the solution (|l^) is stable, while 
the vacuum state is not. In fact, when the initial state of 
the system is a mixture of the two phases: c(x, t = 0) = 
for x < 0, and c(x, t = 0) = c eq for a; > 0, say, then 
the stable phase invades the unstable phase. The front 
between the two phases propagates at a constant speed, 
similar to the case of Fisher waves [0. Here we wish 
to study another inhomogeneous situation, where there 
is a perfectly absorbing trap at the origin instead of the 



initial empty half-space. The trap depletes its immediate 
neighborhood, but a non-trivial steady state is expected 
as the depletion zone created by the trap is continually 
replenished by a stream of particles from the stable phase. 

To derive the appropriate boundary condition, we turn 
back to the discrete representation. In the presence of 
a perfect trap at n = 0, Eq. (||) is then limited to 
< n < m. The special equation for n = is 



dEp t , 
dt 
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Comparison to Eq. (||) yields the discrete boundary con- 
dition 



E-i <m — Eo tm , 
which in the continuum limit becomes 
dE(x,y;t) 



dx 



\x=0 







(13) 



(14) 



In addition, the boundary condition (0) is now replaced 
by 



lim E(0,y;t) = 



(15) 



The stationary solution to Eq. (|j), confined to the 
wedge (0 < x < y), which satisfies the boundary con- 
ditions (§,§), (H) and ©, is 



E s (x, y) = e-rt"-*) + 7(2/ - x)e~' iy 



(16) 



Far away from the trap, as x, y — » 00, this converges to 
the equilibrium result of Eq. ([To|). From (||), we obtain 
the stationary concentration profile: 



c s (z)= 7 (l-e-^) 



(17) 



As expected, there is a depletion zone of size 1/7 = 2D/v 
near the trap, and the concentration grows asymptoti- 
cally to c eq as x — ► 00. 

The IPDF between nearest particles is surprising. 
From Eqs. (||) and ([if]) we obtain the conditional joint 
probability 



P s (x,y) = 7 2 e^(e 7X -l) 



(18) 



Dividing P s (x,y) by c s (x) yields the "forward" IPDF — 
the probability that given a particle at x the next nearest 
particle to its right is at y: 



y-x 



(19) 



The notation chosen here emphasizes the unexpected re- 
sult that p s is translationally invariant. What is more, 
this Poissonian IPDF is characteristic of particles at equi- 
librium, when there are no correlations between their 
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various positions. Indeed, exactly the same IPDF is ob- 
tained for the equilibrium state of Eq. ( |Io| ) , without the 
trap! 

To obtain the "backward" IPDF — the probability 
that given a particle at y the next nearest particle to 
its left is at a distance z — we divide P s (x,y) by c s (y): 
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1y( z ) = 7— 
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(20) 



The fact that q y (z) is not translationally invariant comes 
as no surprise, because of the trap at x = 0. However, 
q y (z) does not normalize properly! The reason for that 
is that there is a finite chance that there are no particles 
between the particle at y and the trap; i.e., the particle 
at y is the nearest particle to the trap. The probability 
that this happens is po(y) = —dE/dy\ x =o, or 



po(y) 



i 2 ye 



(21) 



With this understanding, the proper normalization con- 
dition is 



c s{y) / Qy(z)dz + po(y) = c s (y) 
Jo 



(22) 



which is indeed met. 



single-point density functions. This is achieved by ne- 
glecting all correlations, and by expressing multiple-point 
density functions as products of single-point densities. In 
view of the peculiar IPDF in our problem (Eq. [H]) , this 
seems a promising approximation. 

Let the probability of having a particle at site n 
be p n (t), then the joint probability of having particles 
at both n and m at time t may be approximated as 

Pn,m{t) ~ Pn{t)Pm(t)- In this fashion, our model is de- 
scribed by the equation 



' ' l = ^[-2p n + (1 - p n ){ftn-l + Pn+l)] 

or 



dt 



+ 2^( 1 ~ Pn)(f>n-1 + Pn+l) , 



(25) 



and the boundary condition po(t) — 0. Notice that in the 
stationary limit, and without the trap, Eq. (|25| ) is exact, 
since in the equilibrium state of the infinite system the 
particles really are uncorrelated! Indeed, solving j25| ) 
when dpn/dt = yields c = p/a = v/(2D + av), which 
agrees with c eq of Eq. ( |ll| ) when a — > 0. 

It is tempting to try and connect Eq. (25) to the 



reaction-diffusion approach, by passing to the contin- 
uum limit. This, however, does not work: writing 
c{x,t) = p n (t)/a, and letting a — > while keeping only 
up to first-order terms in a, yields 



IV. REACTION-DIFFUSION EQUATION 

So far, we have presented an exact solution to the 
problem of diffusion-limited reversible coalescence with 
a trap. Exactly solvable models of diffusion-limited reac- 
tions, however, are rare. We now wish to discuss one of 
the most widely used approximation methods, in light of 
the exact results. 

The method is that of reaction-diffusion equations. 
Here one assumes the existence of a mesoscopic length 
scale within which the system is homogeneous and well 
mixed, and where the reaction rates can be accounted 
for as in classical rate equations. At longer length scales, 
variations in the concentration, c(x,t), give rise to dif- 
fusion. In our case, the appropriate reaction-diffusion 
equation is 



dc(x, t) 
dt 



D 



dx 2 



k,c 2 



(23) 



where k\ and k^ represent the "effective" rates of birth 
and coalescence, respectively. The trap at the origin im- 
poses the boundary condition: 



c(0,t) = . 



(24) 



An alternative approach is that of writing down an 
infinite hierarchy of rate equations for the n-point den- 
sity correlation functions, and truncating the hierarchy 
with a Kirkwood ansatz at some convenient stage. In the 
simplest case, one truncates the hierarchy at the level of 



dx 2 
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(26) 



Thus, it is impossible to identify k\ and ki in this man- 
ner, and one is forced to work with the discrete equa- 
tion (p5|). (An approximate solution could be found by 
proceeding with Eq. (26) anyway, ignoring the fact that 
a is supposed to be infmitesimally small.) 

We now use the exact solution of Section III to at- 
tempt and determine k\ and ki. We note that without 
the birth and coalescence reactions the particles would 
simply diffuse with a diffusion constant D — the same 
D as in the hopping rate D/a 2 of the microscopic rules. 

First consider the stationary solution of Eq. (E3|) 
an infinite system without the trap: c eq — kifk^. 



conform with the exact solution of (11), we must have 



ki 
k 2 
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(27) 



The stationary solution to Eq. (g3() with the trap 
the boundary condition (E4J) — is 



^ = |tanh 2 




tanh 



(28) 



where Cqo = c eq is the concentration of particles infinitely 
far away from the trap. 

The concentration profile described by Eq. ( plj ) looks 
similar to the exact result of Eq. (|l7j). One could now 
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use different criteria to further constrain k\ and &2 . De- 
manding the same asymptotic behavior far away from 
the trap; lim^^oo ln[l — c s {x) / c eq ]/x — —v/2D, we get 



_ 3 v * 



(29) 



On the other hand, if we require the same behavior close 
to the trap; (dc s /dx) x= o = (v/2D) 2 , we get 



9f_ 
8 2D 



(30) 



Clearly, it is impossible to fix the short range behavior 
and the long range behavior simultaneously. Instead, one 
may write 



k\ = a ~^~j=j j an d &2 = av 



(31) 



where a is a single fitting parameter of order unity. A 
least square fit in the range < c eq x < 5 is achieved with 
a = 0.889. 



The present model could present an advantage in future 
numerical studies, because of its non-trivial stationary 
state. 

Future work will also include the investigation of two- 
point density correlation functions, as well as the influ- 
ence of a drift away from the trap. Both problems can be 
formulated rigorously within the framework of the IPDF 
method |lC| . Multiple-point density correlation functions 
would shed further light on the breakdown of reaction- 
diffusion equations in low dimensions. Drift away from 
the trap could potentially give rise to an interesting phase 
transition: from a system with a non-trivial steady state 
(such as in our case, when the drift is zero), to a system 
with only a trivial steady state (the vacuum), as the drift 
increases beyond a critical point. 
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V. DISCUSSION 



We have solved the problem of diffusion-limited re- 
versible coalescence with a trap, in one dimension, ex- 
actly. The result is tantalizingly simple: the stationary 
concentration profile is exponential. Moreover, the distri- 
bution of distances between nearest particles (the IPDF) 
is also exponential, similar to that of particles in an equi- 
librium process — the same process in the absence of the 
trap. This does not mean, however, that the distribu- 
tion of particles in the two cases is identical: it is just a 
peculiarity of the IPDF in this particular model. To be 
sure, the distribution of particles in the equilibrium sit- 
uation is fully random and uncorrelated, whereas in the 
presence of the trap it is not! 

We have also contrasted the exact solution with the al- 
ternative, traditional approach of reaction-diffusion equa- 
tions, highlighting the fact that the latter is merely an 
approximation method. Our model provides a clear ex- 
ample where the effective rates of the reaction-diffusion 
equation can be related to the microscopic rates of the 
underlying process, without appealing to renormaliza- 
tion. 

An interesting open question is at what dimension 
reaction-diffusion equations accurately describe the ki- 
netics of the system. Previously, we had conducted nu- 
merical studies of Fisher waves in the coalescence pro- 
cess, suggesting that the critical dimension is d c = 3 flTT| |. 
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